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Introduction 


^he  ancient  Greeks  opened  up  a vast  and  fruitful 
field  for  science  when  their  philosophers  postulated  that 
matter  was  made  up  of  little  indivisible  building  blocks 
which  they  called  "atoms".  Years  passed  and  men  like  I.  -J. 
Thomson, Chadwick,  Anderson,  and  others  showed  that  these 
elementals  particles  were  what  are  now  called  electrons, 
neutrons,  positrons,  etc.  The  word  "atom"  was  preserved 
in  a different  sense  to  mean  the  smallest  particle  of 
matter  that  can  exist.  These  workers  performed  countless 
experiments  to  verify  the  existence  of  these  elemental 
particles.  Early  scientists  in  the  field  of  atomic  physics, 
not  merely  for  the  want  of  anything  better  but  because 
they  really  believed  so,  claimed  that  each  atom  was  a 
miniature  solar  system  with  the  electrons  revolving  about 
the  nucleus  just  as  the  planets  revolve  about  the  sun.  Niels 
Bohr  took  part  of  this  philosophy  over  into  his  theory  of 
the  atom.  He  was  the  first  to  maintain  that  atoms  could 
exist  in  only  certain  definite  energy  states.  His  theory 
embodied  the  principle  that  energy  comes  in  little  discrete 
packets  called  "quanta."  Eventually  his  theory  broke  down 
under  the  weight  of  the  corrections  necessary  to  anply  to 
it  to  make  its  predictions  fit  the  facts,  and  a new  era  of 
scientific  research  came  into  being..  It  was  heralded  in 


II. 


by  the  work  of  Schroedinger  and  others  who  founded  the 
branch  of  mathematical  physics  known  as  wave  mechanics. 

The  application  of  wave  mechanics  to  the  practical 
world  of  solids  which  we  can  see  and  touch  rather  than  to 
the  intangible  world  of  electrons,  protons,  etc.  seemed  to 
the  writer  to  be  a field  which  should  prove  tremendously 
interesting.  Accordingly  when  severe  doubt  was  cast  on 
the  method  employed  by  Fuchs  to  calculate  the  cohesive 
energy  and  lattice  constant  of  Cu,  it  seemed  an  excellent 
opportunity  to  enter  the  field.  Consultation  with  others 
confirmed  this  feeling  and  it  was  decided  to  test  Fuchs* 
method  by  applying  it  to  a considerably  heavier  element. 

A.g  was  chosen  for  reasons  explained  in  more  detail  in  the 
body  of  this  paper.  The  work  proved  difficult  and  extremely 
tedious  as  methods  of  successive  approximation  are  apt  to 
be.  Furthermore,  the  necessity  of  including  the  exchange 
energy  between  the  valence  electron  and  the  core  electrons 
at  least  quadrupled  the  amount  of  necessary  calculation. 
Fortunately  as  the  different  problems  connected  with  the 
various  phases  of  the  solution  presented  themselves  means 
were  found  to  solve  them.  The  actual  final  results  were 
highly  satisfactory  and  on  the  average  more  accurate  than 
those  of  Fuchs,  thus  completely  answering  the  purpose  of 
the  dissertation. 

The  writer  would  like  to  thank  Professor  Foysl  M.  Frye 
for  his  valuable  suggestions  and  inspiration  throughout  the 


III. 


solution  of  the  problem,  Professor  J.  C.  Slater  of  the 
Massachusetts  Institute  of  Technology  for  his  advice  on 
methods  of  solving  the  differential  equations  involved, 
and  Miss  Mary  F.  Wood  for  help  in  preparation  of  the 
manuscript . 
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Part  I 


Review  of  the  Field  and  Methods  Employed 


The  basis  of  this  investigation  is  the  wave  mechanics 
as  developed  by  Schroedinger^in  1926  and  now  embodied  in 
his  justly  famous  wave  equation  as  given  below: 


The  form  of  this  equation  most  valuable  to  us,  however,  is 
the  differential  equation  for  the  radical  function 


This  paper  deals  particularly  with  the  application  of 

this  equation  to  solids  and  the  subsequent  determination 

by  it  of  the  cohesive  energy  and  interatomic  distance  of 

silver.  The  Schroedinger  Equation  (henceforth  abbreviated 

"S  - equation”)  has  been  applied  successfully  to  isolated 

atoms  for  the  past  fifteen  years  but  its  effective  application 

2 

to  solids  dates  from  the  work  of  Wigner  and  Seitz  in  1933. 
Hence  for  purposes  of  reviewing  the  field  to  date,  it  will 
be  logical  to  use  their  work  as  a starting  point. 

Wigner  and  Seitz  chose  for  their  initial  task  metallic 
Na.  This  was  not  a chance  choice  from  the  list  of  metals 


1«  E.  Schroedinger,  Ann.  D.  Phys.,  79,  361,  487;  80, 
437,  81,  109  (1926) 
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2.  E.  Wigner  and  F.  Seitz,  Phys . Rev . , 43,  804  (1933) 
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but  a careful  one  based  on  the  following  considerations: 

a.  Prokofjew  had  already  worked  out  a series  of 
analytic  functions  for  the  potentials  in  Na. 

b.  It  had  comparatively  few  electrons. 

c.  It  had  a large  interatomic  distance,  hence  more 
closely  approximating  a set  of  isolated  atoms  than  say 
the  noble  metals. 

The  method  of  solution  embodied  what  is  now  known  as 
the  "cellular  method."  This  method  consists  in  constructing 
about  each  atom  in  the  crystal  a polyhedron  formed  by  planes 
bisecting  perpendicularly  the  lines  joining  a particular 
atom  with  its  nearest  and  in  some  cases  next  nearest  neighbors. 
Figure  1.  illustrates  the  method  of  division  in  two  dimensions. 


! 


i 

Figure  1 


The  purpose  of  the  perpendicular  bisection  was  that  at  the 


3.  W.  Prokofjew,  Zeits  f.  Physik. , 58,  255  (1929) 
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mid-points  the  electric  forces  from  each  atom  would  cancel, 
hence  since  the  boundary  surfaces  will  he  neutral;  we  can 
more  or  less  ignore  all  save  one  polyhedron.  This  concentrating 
on  one  polyhedron  alone  is  quite  effective  in  metals  with 
large  interatomic  distance  as  in  the  alkali  metals  but  in 
the  noble  metals  it  cannot  be  done  safely.  It  has  been 
found  that  if  we  replace  the  polyhedra  by  spheres  of  equal 
volumes  that  the  error  incurred  will  be  negligible  compared 
with  subsequent  errors  to  be  discussed  later. 

In  an  isolated  atom  we  know  that  a function  is  not 
a repeating  one  but  Bloch^has  shown  that  in  a solid  the  ^ 
is  a repeating  function  of  the  form 


% e UK 


Thus  the  curve  must  be  smooth  as  it  passes  from  one 
sphere  to  another.  If  the  radius  of  the  sphere  of  maximum 
stability  is  r , then  Bloch* s boundary  condition  takes  on 
the  form 

- /£L'° 


It  has  proved  much  more  convenient  to  use  the  S - equation 
as  given  in  (2)  rather  than  in  (l)  hence  (4)  becomes 


4.  B.  Bloch,  Zeits  f.  Physik,  52,  555  (1928) 
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This  equation  when  differentiated  and  solved  for  dR  gives 

\dH  % I ar 

6-  U*  ' * 

In  other  words,  the  value  of  r at  which  the  slope  of  the 

R curve  is  R is  equal  r . This  can  be  found  easily  by 
r 0 

drawing  tangents  from  the  origin  to  the  R curve  at  distances 
from  the  nucleus  greater  than  that  of  the  outermost  bound 
shell.  In  Na  this  would  be  the  2p  shell,  in  Cu  the  3d 
shell,  and  in  Ag  the  4d  shell.  It  will  be  noticed  that 
for  R curves  corresponding  to  energies  numerically  less 
than  the  energy  of  the  sphere  there  will  be  two  points  of 
tangency  and  at  energies  greater  no  points  of  tangency. 

The  presence  of  the  two  points  is  easily  explained  as  they 
correspond  to  the  energy  the  sphere  would  have  if  compressed 
or  expanded.  Since  neither  of  these  conditions  is  a stable 
one,  that  value  of  the  energy  which  gives  only  one  point  of 
tangency  is  the  energy  of  the  sphere  in  its  most  stable 
position.  This  also  shows  up  as  the  minimum  of  the  energy 
vs.  rQ  curve. 

Having  now  in  mind  the  general  plan  and  boundary 
conditions,  it  remains  to  solve  the  S - equation  for  the 
valence  electron  since  it  was  assumed  that  the  core 
electrons  are  bound  into  rigid  shells.  This  assumption 
is  necessary  as  it  would  be  practically  impossible  to 
treat  all  of  the  electrons  in  the  manner  to  be  outlined. 


. 


•• 
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Spherical  symmetry  is  assumed  for  the  valence  electron. 

Hence  we  can  use  the  S - equation  of  the  form  in  equation  (2). 

The  first  step  in  the  solution  of  an  S - equation  is  the 

establishment  of  the  value  of  V (r  ),  the  pseudo-classical 

potential  energy.  One  must  say  "pseudo"  since  it  is  actually 

the  potential  energy  field  in  which  the  electron  finds 

itself  and  since  the  potential  is  a combination  of  effects 

due  to  the  nucleus  and  core  electrons,  it  cannot  in  most 

cases  be  successfully  approximated  by  a Coulomb  potential 

term.  Wigner  and  Seitz  used  a half  empirical  set  of 

5 

potentials  derived  by  Prokofjew.  In  addition,  they 
realized  that  it  might  be  possible  for  more  than  one  valence 
electron  to  be  inside  the  same  sphere  at  the  same  time. 

This  they  assumed  could  not  happen  and  hence  took  into 
consideration  in  part  at  least  the  correlation  between 
electrons.  Having  these  Prokofjew  potentials,  the 
S - equation  was  solved  by  the  method  of  finite  differences 
for  different  values  of  W and  the  values  of  R corresponding 
to  each  plotted.  Curves  similar  to  that  shown  in  Figure  2 
were  obtained.  From  these  as  already  described  the  values 
of  rQ  were  gained  and  the  W vs.  rQ  curve  plotted.  Up  to 
now,  however,  the  fact  that  the  valence  electrons  were  in 
motion  v/as  neglected,  but  now  the  energy  due  to  motion 
must  be  taken  into  consideration  in  order  to  get  the 
total  energy  of  the  valence  electron.  This  was  necessary 


5.  W.  Prokof jew,  Zeits  f.  Phys.,  58,  255  (1929) 
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since  the  cohesive  energy  is  defined  as  the  ionization 
energy  minus  the  total  energy.  If  one  assumes  that  the 
velocities  of  the  electrons  are  distributed  according  to 
Fermi -Dirac  statistics  and  that  the  electrons  are  practically 
free,  then  this  additional  energy  can  be  calculated  from 
the  formula 


The  ionization  potential  may  be  also  obtained  analytically 
by  observing  that  the  energy  curve  approaches  an  asymptotic 
value  as  rQ  increases.  This  is  easily  interpreted  as  the 
energy  required  to  expand  the  valence  electron  (the  core 
electrons  are  rigid)  to  infinity  or  in  other  words,  the 
ionization  potential.  The  energy  curves  for  Na  are  given 
in  Figure  2. 
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Figure  2 
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In  the  course  of  their  investigations  Wigner  and  Seitz 

O 

together  and  Wigner  alone  attempted  to  take  into  account 


7.  E.  Wigner  and  F.  Seitz,  Phys . Rev . , 46,  509  (1954) 
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E.  Wigner,  Phys . Rev . , 46,  1002  (1934) 
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the  facts  they  knew  had  been  omitted  in  their  first  treatment 
based  on  Prokofjew  potentials.  The  corrections  they  made 
were  based  on  the  Hartree  field  for  Na  and  are  given  below: 

a.  Interaction  of  electron  with  electron  in  its  own 

sphere  with  no  correlation  6e^ 

ro 

b.  Decrease  in  energy  due  to  the  "Fermi”  hole 

between  electrons  with  parallel  spins  458e^/r0 

c.  Decrease  due  to  hole  for  electrons  with  anti- 

paralled  spin  - be2/r0 

b is  a function  of  rQ 

d.  Van  der  Waal's  attraction  between  ions  .0024  e.v.  /atom 

e.  The  repulsion  between  shells  .Ole.  v.  /atom 

In  Table  I the  results  of  Wigner  and  Seitz  calculations 
before  and  after  the  several  corrections  had  been  applied 
are  compared  with  the  experimental  data 


Table  I. 

Binding  Energy 

(in  Kilo-cal/gm  atom) 

Lattice  constant 
(in  Angstroms) 

Comnressibility 
(in  Kg-cm^x  10°) 


Before 

Correction 

23T5 


4.20 


16 


After 

Correction  Experimental 

£BTI 25T5 


4.62  4.25 

16  10 


This  method  then  can  be  said  to  be  successful  but  the 


number  of  corrections  necessary  gave  an  inkling  that  the 


9.  D.  R.  Hartree,  Proc.  Camb.  Phil.  Soc.,  24,  111  (1928) 
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Prokof jew  potentials  were  probably  not  complete.  This 
suspicion  was  later  verified  by  Gorin^/ho  attempted  to  construct 
a set  of  potential  functions  for  K by  Prokof jew's  method. 

He  found  that  it  was  impossible  to  do  so  at  all  accurately. 

The  reason  for  the  failure  was  felt  to  lie  in  the  possibility 
that  the  exchange  terms  could  not  be  replaced  by  ordinary 
Coulomb  potential  terms  to  even  a moderate  approximation. 

Realizing  that  the  Prokof jew  method  was  not  yielding 

good  results,  Gorin  turned  to  the  Hartree  ion-core  field 

for  K which  had  already  been  worked  outl1  Here  again  he  was 

unfortunate  in  that  his  direct  results  came  out  very  poorly, 

the  cohesive  energy  coming  out  6 kilo-cal/mol  instead  of 

22-6  kilo-cal/mol.  This  time  the  error  was  claimed  to  be 

in  the  Hartree  field.  The  error  in  the  lowest  level  of 

the  atom  being  .735  e.v.  without  exchange  and  .347  e.v.  with 

exchange.  In  the  cases  of  Na  and  Li  the  Hartree  field  had 

worked  rather  well  so  that  it  is  generally  now  believed 

that  although  the  Plartree  functions  are  admittedly  inaccurate, 

they  are  the  best  starting  functions  available.  Gorin  did 

try  to  "D.octor  judiciously"  his  results  by  inserting  a 

factor  . 388  -t'. 347 , the  .338  being  the  exchange  energy 
.338 

between  the  valence  electron  and  the  core  electrons  and 


10. 


11. 


E.  Gorin,  Physic  Z.  Saw j , 9,  328  (1936) 

D.  R.  Hartree,  Proc.  Roy.  Soc.,  143,  506  (1934) 


the  .547  is  the  corresponding  correlation  energy.  Table  II 
gives  a summary  of  his  results. 

Table  II 
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Cohesive  Energy  in  kilo-cal/mol 


First  Calculation 

6.0 

Corrected  Calculation 

14.3 

Observed 

22.6 

Both  Wigner  and  Seitz  and  Gorin  hiad  worked  with  alkali 

12 

metals  and  it  remained  for  Fuchs  to  attempt  the  first  of 

the  noble  metals.  These  metals  are  characterized  by  cores 

which  are  much  larger  than  those  of  the  alkali  metals  in 

comparison  with  their  respective  atomic  volumes.  Thus 

the  possibility  of  a much  stronger  ion-ion  interaction  was 

suggested.  In  addition,  the  noble  metals  exist  as  face 

centered  cubic  crystals  whereas  the  alkalis  are  body  centered. 

Fuchs  chose  Cu  from  the  noble  metals  since  it  was  the 

simplest  of  the  group  and  the  wave  functions  for  Cu  had 

13 

been  previously  worked  out  by  Hartree.  The  cellular  method 
developed  by  Wigner  and  Seitz  was  used  with  the  exception 
that  a Hartree  self-consistent  field  in  place  of  the 
semi-empirical  Prokofjew  field,  This  exception  necessitated 
using  the  Fock  form  of  the  radial  S - equation  for  the 
valence  electron  since  the  Hartree  functions  do  not  contain 


12.  K.  Fuchs,  Proc.  Roy.  Soc.,  151,  585  (1935) 

13.  D.  R.  Hartree,  Proc  Roy.  Soc.  A,  141,  282  (1933) 


143,  506  (1933) 
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the  exchange  energy  implicitly  as  do  the  Prokofjew  functions. 
This  equation  for  the  valence  electron  (the  4S  electron  for 
Cu)  is  of  the  form. 


and  the  summation  is  over  all  the  inner  electrons  with  the 
same  spin  as  the  valence  electron  (i.e.  half  the  number  in 
each  shell).  TheJ^  S are  the  function  for  the  closed  shell 
electrons  and  are  given  by 


and  the  spherical  harmonic  part  is  normalized  as  shown  in 
equation  11. 


11. 


Q J & </  0 - 


Y rr 


The  valence  electron  function  is 


In  his  work  Fuchs  showed  that  G(r,rf)  may  be  put  in 
the  following  form  which  is  made  more  convenient  for  numerical 
calculation 


15. 


- 


. 


. 

. - 
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and  n and  1 are  the  principal  and  azimuthal  quantum  numbers 
respectively. 

As  a result  of  using  the  Fock  equation,  the  calculations 
were  much  more  difficult  and  tedious  than  those  based  on 
the  simuler  equation  used  by  Wigner  and  Seitz.  The  method 
of  solution  adopted  was  to  neglect  the  exchange  term  in  (8) 
and  solve  the  remainder.  This  gave  a fairly  good  value  for 
f(r).  Then  the  exchange  integral  was  calculated  using 
Hartree's  functions  for  the  core  electrons  and  the  approximate 
value  of  the  4S  function  just  mentioned.  The  result  was 
then  substituted  into  equation  (8)  and  the  subsequent 
inhomogeneous  equation  solved.  Fuchs  felt  that  the  result 
of  one  approximation  obtained  in  this  manner  was  of  sufficient 
accuracy  not  to  warrant  recalculation.  The  Fermi  energy  was 
obtained  from  the  formula 


14 


where  U - atomic  volume  - 47Tr03 

g 

His  first  results  are  given  in  Table  III 


Table  III 


Lattice 


Lattice  Heat  of  Vaporization  Compr^ssibi 

Constant  in  A in  Kilo-cal/mol Kg/cmLx  10~ 


Theoretical 


4.1 


•41.4 

3 

71.0 


2.6 


Experimental 


3.6 


70 


These  results  were  obtained  in  precisely  the  same  manner  as 
those  of  Wigner  and  Seitz  which  have  already  been  explained. 


12. 


Since  the  compressibility  was  very  much  in  error,  Fuchs 
calculated  in  addition  the  ion-ion  repulsion  and  found 
that  it  was  as  had  been  suspected  much  larger  in  the  noble 
metals  than  in  the  alkalis.  An  additional  correction  was 
attempted  as  regards  the  Fermi  energy  as  determined  by  (14). 
That  equation  was  derived  on  the  basis  of  perfectly  free 
electrons  and  the  valence  electrons  in  the  lattice  cannot 
be  said  to  be  "perfectly  free."  The  expression  was  corrected 
to 


and  g is  the  vector  drawn  from  the  origin  in  K-space  to 
the  nearest  point  of  the  first  zone  boundary  for  the  face 
centered  cubic.  Fuchs]-4using  approximation  methods,  showed 
that  the  value  of  obtained  from  (15)  was  very  close  to 
that  given  by  equation  (14).  Seitz  has  since  pointed  out 
that  the  approximations  Fuchs  made  in  this  last  instance 
were  unjustified.  Fuchs  results  after  the  corrections  had 
been  applied  are  given  in  Table  IV. 


14.  K.  Fuchs,  loc.  cit,  (12)  p.  509 


15. 


F.  Seitz. Modern  Theory  of  Solids,  p.  570 
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Table  IV 

Lattice  Heat  of  Vaporization  Compressibility 

Constant  in  A in  Kilo-cal/mol  Kg^/cm^x  10° 

Theoretical  4.2  34.1  .69 

Experimental  3.6  71.0  .70 

As  regards  Fuchs  work,  we  may  say  that  his  results  while 
not  as  accurate  as  those  of  Wigner  and  Seitz  are  of  consider- 
able importance  in  that  it  served  as  a preliminary  comment 
on  Hartree* s functions  when  used  in  calculating  the  constants 
of  the  noble  metals.  Fuchs  himself  remarked  on  this  point. 

When  he  pointed  out  that  4 % of  the  charge  distribution  of  the 
3d  shell  for  Cu  lay  outside  the  neutral  polyhedron.  The 
continuance  of  this  investigation  on  the  use  of  Hartree 
functions  for  even  heavier  metals  forms  the  basis  of  the 
original  part  of  this  paper. 

Before  we  discuss  the  thesis  proper,  it  will  be 
particularly  advantageous  in  the  light  of  the  development 
to  date  and  in  anticipation  of  that  which  is  to  come  to 
discuss  in  more  detail  the  Hartree  and  Fock  method  for 
solving  the  S-  equation.  First  let  us  assume  that  the 
for  the  whole  atom  can  be  represented  as  a product  of  one 
electron  functions  such  as 


Hartree  suggesteJ6that  each  one  electron  function  in 
(17)  should  satisfy  a one  electron  S - equation  where  the 
V must  take  into  consideration  not  only  the  effect  of  the 


16.  D.  E.  Hartree,  Proc . Cambridge  Phil.  Soc.,  24,  105  (1926) 
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charge  distribution  of  the  nucleus  but  also  that  of  all  of 
the  other  electrons  save  the  one  being  considered  as  well. 


summation  term  arising  from  the  other  electrons,  and  £.  is 
the  energy  parameter.  The  summation  term  takes  care 
partially  of  the  electron-electron  interaction  and  the 
prime  on  it  indicates  that  i cannot  = j . This  takes  care 
to  a certain  extent  of  the  Pauli  Principle  but  incompletely 


of  the  position  of  the  other  electrons.  This  causes  a 
definite  error  to  arise  since  obviously  cannot  be 
independent  of  the  position  of  the  other  electrons.  Thus 
if  we  consider  He,  if  one  electron  is  on  the  right-hand 
side  of  the  nucleus,  the  other  electron  would  more  probably 
be  on  the  left-hand  side  than  the  right  due  to  the 
electrostatic  repulsion  of  like  charges.  This  effect  is 
disregarded  completely  in  the  Hartree  method. 


Hartree1 s method  for  determining  the  T is  one  of 
successive  approximation  and  is  known  as  "the  method  of  the 
self-consistent  field."  The  name  arises  from  the  type  of 
approximations  used.  To  illustrate  very  briefly  the  method, 
Hartree  uses  in  solving  an  atom  let  us  consider  He»  The 
S - equations  for  the  two  electrons  of  He  are 


The  equation  for  each  electron  then  becomes 


17 


18. 


where  V.  is  the  field  arising  from  the  nucleus,  the 
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In  solving  these  equations  one  estimates  and  then 
solves  for  ^ by  numerical  integration.  Then  Vg  may  be 


approximated  by  using  the  formula 


e 


81  • ^ 


*-  © 


Now  knowing  Vg,  r ^ may  be  found  and  a new  approximation 

made  to  V-,  ( r ) . This  is  kept  up  until  Vn(r) -Vn-l(r)  for 

Vn(r) 

V is  less  than  any  given  ratio  which  point  the  assumed 

and  calculated  fields  are  consistent  with  each  other. 

Foci3realized  that  Hartree’s  one  electron  approximation 

was  accurate  only  if  the  electrons  were  all  in  the  same 

19 

state  or  all  free  electrons.  Thus  he  and  Slater  invented 
a set  of  equations  based  on  the  variation  principle  rather 
than  on  a Hartree  principle.  Briefly  their  theory  stated 
that  a better  method  for  obtaining  the  best  for  the 
atom  was  to  use  those  which  made  the  integral 

/r 
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A minimum  for  any  small  variation  in  ^ 
the  normalizing  restriction  that 


consistent  with 
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H is  the  Hamilton  operator  and  is  given  by 


24 


. H < H;  * j £ i!  ( i. j - / • • -~) 


'■j 


where  depends  only  upon  the  variables  r^  and  is  the 

same  function  of  these  as  H.  is  of  r-.  The  result  of 

3 ^ 

employing  the  variation  principle  is  to  change  the  form 
of  the  y function  for  the  atom  from  that  given  in 
equation  (17)  to 

$ (*,)  0,  -\  ) ■ 4>,  ’ 

v ’ -h  


25. 





In  which  the  0 are  the  functions  of  the  electrons.  The 

1 c 

general  form  of  Fock* sSOequation  using  this  type  of  / i; 


- 

26  . eL~>n, 


J 


Ac 


The  /\  cj  are  the  parameters  necessitated  by  imposing  the 
orthoganility  condition 

jfrfydT 4,  y,J,  y ■)  . oJTy 

where  ^represents  the  apin  coordinate. 

In  (19)  and  (20)  the  Hartree  equations  were  written  but 
for  H and  it  is  interesting  to  compare  them  with  the  Fock 


20.  F.  Seitz,  Modern  Theory  of  Solids,  p,  244 
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equations  given  below  for  the  same  element.  (Note  the 
increased  complexity) 


Equations  (27)  and  (28)  were  written  in  atomic  units. 
The  reason  that  Fock’s  equations  are  superior  to 


those  of  Hartree  is  that  the  determinental  form  for  which 


Fock  uses  allows  for  the  possible  antisymmetry  in  the 

coordinates  of  the  different  electrons.  It  is  interesting 

to  note,  however,  that  Focks  equations  can  be  reduced 

exactly  to  Hartree Ts  equations  if  one  assumes  the  one 

electron  approximation  as  did  Hartree.  Because  of  the 

extreme  complexity  of  Fock's  equations  comparatively  few 

21 

of  them  have  been  solved  and  Hartree  functions  are  usually 

2? 

employed  particularly  in  work  on  solids.  In  addition  Fock 


21.  Nfi,  Lj_,  V.  Fock  and  Petrashen,  Phys.  Z.  Sow.jet, 
6,  368  (1934);  8,  547  (1935)  K,  D.  P.  Hartree,  Proc.  Poy. 
Soc . , 166,  450,  (1938)  Cu,  P.  P.  Hartree,  Proc.  Poy.  Soc. 


157,  490  (1936) 

22.  V.  Fock,  Zeits  f.  Phys.,  81,  195  (1933) 
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has  shown  that  for  a valence  electron  the  radial  equation 


can  be  put  in  the  form 

-JL  dX 


y- 


JvJ  - 


/v  c O 


This  equation  is  much  the  same  as  Hartree’ s writh  the 
exception  of  the/  This  integral  is  called 


the  "exchange  integral"  and  enters  because  the  ^ for  the 
atom  was  a determinant.  Equation  (29)  can  be  solved  by- 
using  Hartree  functions  for  the  core  electrons  as  will 
be  shown  in  the  thesis  proper. 

A detailed  derivation  of  Hartree fs  and  Eock's 
equations  is  not  necessary  for  the  purposes  of  this  paper 
but  may  be  found  in  standard  texts  on  this  subject. 


25.  E.  Seitz,  Modern  Theory  of  Solids,  pp.  670  - 680  inc. 
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Part  II 


Calculation  of  the  Cohesive  Energy  and  Lattice  Constant 
of  Silver  by  the  Method  of  Fuchs. 

A.  Purpose: 

In  view  of  the  fact  that  Fuchs  and  Wigner  and  Seitz 

had  to  apply  several  correction  terms  to  their  first 

determinations  from  the  S - equation  in  order  to  bring 

their  results  into  closer  agreement  with  the  experimental 

values,  it  was  felt  that  an  extension  of  Fuchs’  method 

to  an  even  heavier  element  might  experience  a complete 

breakdown.  The  breakdown  referred  to  would  be  evident  by 

results  of  an  entirely  wrong  order  of  magnitude.  Accordingly 

Ag  was  chosen  since  it  is  53.9$  heavier,  has  47  electrons 

as  compared  with  the  29  of  Cu,  and  Hartree  functions  and 

24 

potentials  were  available. 

Thus  the  purpose  of  this  dissertation  is  to  test  the 
possibility  of  applying  Fuchs’  method  to  solids  heavier 
than  Cu. 

B . Method: 

As  this  thesis  is  a test  of  Fuchs’  method,  it  wrill  be 
employed  of  course  explicitly  in  determining  the  cohesive 
energy  and  lattice  constant  of  Ag.  The  method  was  outlined 


24.  M.  M.  Block,  Mem,  and  Proc . Manchester  Lit.  Phil. 
Soc . , vol.  59,  no.  3,  28  (1934-35) 
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in  some  detail  in  Part  I of  this  paper.  Since  this  method 
differs  from  the  more  well  known  method  of  Wigner  and  Seitz 
mainly  in  the  use  of  Hartree  functions,  it  will  he  of 
interest  to  note  .just  how  these  functions  alter  the 
numerical  calculations  involved.  First  of  all,  as  has 
been  mentioned  before ,sthe  use  of  Hartree  functions 
necessitates  the  introduction  of  an  integral  called  the 
"exchange  integral"  into  the  S - equation.  The  equation 
Wigner  and  Seitz  used  and  the  one  Fuchs  used  are  given  in 
equations  (30)  and  (31)  respectively  for  purposes  of 
comparison. 

30.  £±.(±4-0$  '° 

CL  Sv  v 

si.  tJjx (fc-6)<l>+  x f*  q ,0,  J <Pw)Ja,  - o 

d a>  v - o 

These  equations  where  G(rr')  - Xi ^ ^ <r ' 

Jo  c JAj  -sJj 

represent  considerable  simplification  as  can  be  seen  by 
writing  equation  (50)  in  the  more  customary  form 


32. 
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25.  cf.  pp.  9-13 

26.  cf.  p.  10 
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The  simplification  was  effected  by  using  atomic  units  given 
in  Table  V. 

Table  V 

Unit  of  length  = h2/4  7T  me^ 

Unit  of  charge  - c = 4.77  x 104^  e.s.u. 


Unit  of  mass  - me  - 9.035  x 10 


-28 


gms 


Unit  of  energy  = 2 7rme4  = 13.53  e.v. 

h^ 

The  W in  equation  (32)  is  thus  measured  in  Rydbergs  and 
for  convenience  sake  will  be  set  equal  to  -g  v is  set 

^ y 

equal  to  -Vc . Nov/  making  these  substitutions  and  multiplying 
the  equation  (32)  through  by  a factor  2,  one  obtains 


33. 


TT^-th.  -e?  ol  Su 


u)<t 
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If  the  operations  indicated  are  performed  and  the 
atomic  units  used,  equation  (35)  becomes  equation  (30).  It 
is  necessary  to  note  again  that  the  EU  in  equation  (30),  (31) 
and  (53)  is  in  double  Rydberg  units,  twice  the  ionization 
potential  of  hydrogen.  The  ensuing  calculations  were 
considerably  eased  by  this  scheme. 

The  method  of  solution  adopted  was  to  ignore  the 
exchange  integral  in  equation  (31)  and  solve  the  remainder. 
This  gave  an  approximate  value  of  <f)  , the  function  for  the 
valence  electron  (the  5S  electron  in  the  case  of  Ag) . 


27.  Mott  and  Jones,  Properties  of  Metals  and  Alloys,  p.  519 
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This  approximate  (p  for  an  alkali  metal  is  quite  close  to 
the  correct  one  viiereas  in  the  noble  metals  the  (j)  obtained 
in  this  manner  and  that  obtained  by  solving  the  complete 
equation  are  very  different.  This  is  shown  in  Figure  3 
where  curves  for  £=  2.1  with  and  without  exchange  for  Ag 
are  plotted. 

The  next  step  was  to  calculate  the  value  of  the 
exchange  integral  and  it  is  this  tern  which  increases  the 

OO 

difficulty  of  solution  immensely.  Fuchs  has  shown  that 
the  term  G-(r,rf)  may  be  put  into  the  form 


34.  (j(Sb  a') 


K 


1 


In  this  expression  the  f(r)  and  f(rf)  are  the  functions  of 
the  core  electrons  given  by  Block  and  normalized  to  . 

It  must  be  noted  that  the  summation  extends  over  all  core 
electrons  having  the  same  spin  as  the  valence  electron. 
Thus  the  exchange  integral  can  be  put  into  the  form 


J 


35. 

36.  If  we  set  // 

37.  And  then 


x.i  'Vv 


u 


Aj 


A,  A 


.A,  A 


At 


The  exchange  integral  may  now  be  put  into  the  form 


28.  K.  Fuchs,  Proc.  Boy.  Soc..  151,  587,(1935)  * 
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This  integral  can  be  solved  by  graphical  integration.  The 

value  of  the  integral  at  a particular  r is  of  prime  importance 

29 

because  Block's  table  of  the  term  2VC  which  appears  in 
equation  (31)  is  also  given  not  as  an  analytic  function  of 
r but  as  in  Table  VI. 

Table  VI 


r 

P-! 

tSJ 

C\2 

2VC 

1.4 

5.09 

3.64 

1.6 

4.04 

2.53 

1.8 

3.38 

1.88 

2.0 

2.94 

1.47 

The  2ZT)listed  in  Block's  tables  is  equal  to  r times  2VC 
where  V is  the  potential  for  Ag.  Then  in  solving  the 
S - equation  the  value  of  2VC  at  a particular  r is  taken 
and  the  corresponding  value  of  the  exchange  integral 
applied  to  complete  the  equation.  For  example  the  S - equation 
for  the  valence  electron  at  a distance  r=1.0  and  an  energy 
equal  to  2.1  double  Rydbergs  is 

39.  CA  + (lit.  - J.A  <j>  + Ax.°oo^oV  = ° 

a /i/ 

To  obtain  the  value  of  the  exchange  integral  at  a particular 


29.  loc.  cit. 
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r one  must  first  expand  A'r,rf  setting  r=1.0  or  A'-j  Q,rf. 
The  actual  expansion  of  A'-^Qjr*  is  given  in  Table  VII. 

To  appreciate  the  work  that  is  involved  in  calculating  the 
exchange  integral  one  must  realize  that  fifteen  such 
expansions  were  performed  and  this  is  only  the  first  step. 
In  Table  VII  the  functions  in  the  expansion  of  A1 ^ q>.2 
and  A* ,.4  are  indicated  symbolically  to  show  the  scheme; 
the  other  values  of  Af]_.o,rf  are  expanded  and  the 
numerical  values  of  the  functions  inserted  directly. 


Table  VII 

A,1.0,.2"1X-2  f3,0(1-0)f3,0(-2)+3x-2  f3,l(1-°)f3,l(-S)* 
5x.2f  (1.0)f  (.2)i-lx.2  f (1.0)  f4  0(.2)« 

3x.2f4  (1.0)  f4  (.2)+5x.2f4  (1.0)  f4  g(.2) 

A’l.0,.4--lx-4  f3,0(1-0)  f3,0(-4)T3x-4  f3,l(1-0)  f3,l 
(.4)i-5x.4  ts  ,(1.0)  f3  g(.2)+lx.4  0(1.0) 

f4  0(.4)-t3x.4  f4  1(1.0)  f4  1(.4)-t5i.4  f4  g 

(1.0)  f (.4) 

A*  =1x1/. 6(  .096107)  ( .S3244)+3x(  .10307)  ( .84745)*- 

1.0, .6 

5x. 6 ( .12559)  ( .76423 ) t-lxl/6  (1.1182)  ( .73754)  + 
3x( 1 . 0918) ( .60710) r5x. 6( .83904) ( .39166) - 
-lxl/8(  .096107)  ( . 31453 )-*■  3x(  .10307)  ( .32066)-* 

5x. 6 ( .12559) ( .32600  ) +1x1/8(1.1182) (1.11896) * 
3x( 1 . 0918) ( 1 . 0660 ) i 5x. 6 ( .83904) ( .71859) 


1.0, .8 


1.0, 1.0 


('lx(  .098107 )2t  3x(  . 10307 )2t  5x(  . 12559  )2t  lx(  1.1182) 
3x(l. 09184) 5x( .33904)2 


* 


. 


. 


« 





m 
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At1.0  i.2~lx(  .096107)  ( .026221)  + 1x3/1. 2V(  .10307)  ( .034356)+ 

1. 0x3/1. 2*t  .12559)  ( .045426 ) +lx( 1.1182 ) ( .36443)  i- 
lx. 3/1. 2(1. 0918)  ( .92713)  *1x5/1. 2*(  .33904)  ( .33904) 

A»  n,  ,rlx( .096107) ( .0087371)^1x3/1. 4( .10307X .011452) * 

1x5/1. A .12559) ( . 016033) * lx ( 1 . 1132 ) ( .60470)+ 
1x3/1.4(1.0913) ( .72241) +1x5/1. 4t .85904) ( .78496) 

A,1.0  1.6-  0 + 0 , 1x5/1. 6^  .12559) ( . 0053443)+ lx( 1 . 1182 ) ( .40644)* 
1x3/1.6(1.0918)  ( .54122)  + 1x5/1. 6l[  .33904)  ( .70957) 

A?1.0  1#3T  1*5/1. 8*(  .12559)  ( .002.672)-rlx(l. 1182)  ( .26567)* 

1x3/1.8(1.0918) ( .38591)+ 1x5/1. 8^ .83904) ( .63010) 

At1.0,2.0-lx(1*1182) ( .16852)+ 1x3/2. 0(1. 0918) ( .27061) * 

1x5/2. 0\  .33904) ( .55589) 

A'l  .0  g.g"lx(  1-H82)  ( .10706)+ 1x3/2. 2(1. 0918)  ( .13825) + 

1x5/2. 2^ .83904) ( .48425) 

A * 1 . 0 , 2 . 4~ lx ( 1 * 118^ ) ( .065427)  + lx3/2.4  (1.0918) ( .12707)  + 

1x5/2. 4*[  .83904)  ( .42034) 

A f 1 . 0 , 2 . 6 " lx ^ 1 * 118^ ) ( .039653) *1x3/2. 6 (1.0918  ) ( .087065 ) + 

1x5/2. .83904) ( .36462) 

A* 1 . 0 , 2 .3" lx( 1 • 1182 ) ( .025774) +1x3/2. 3(1. 0918) ( .058828)+ 
1x5/2.3^.83904) ( .31464) 

A*  j ,0>3 . 0 ”lx ( 1 • H82 ) ( .015861)+  1x3/3.0(1.0913)  ( .040003)  + 

1x5/3. o\ .83904) ( .27121) 


-r-  n ... 
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Table  VIII  is  the  complete  set  of  values  of  A,rrt  calculated 
for  r varying  from  r=0.0  to  r-3.0  by  steps  of  Ar  = . 2 ■ 

Table  VIII 


r’ 

A* .2,r* 

A’ . 4 , r ’ 

At .6,r* 

A f ^ . 

K .8,rT 

A'l.0r* 

A f 1 . 2 , r 

tAt1.4}r 

. 2 

124.41 

33.464 

6.4435 

1.3136 

2.7198 

2.3482 

1.6878 

.4 

33.464 

56.169 

15.258 

2.3545 

.848786 

1.2295 

.97768 

.6 

6.4435 

15.258 

13.784 

7.8587 

5.0402 

3.3379 

2.2061 

.8 

1.3136 

2.3546 

7.8575 

10.431 

7.8668 

5.5999 

5.7225 

1.0 

2.7198 

.84876 

5.0402 

7.8668 

8.4666 

5.9728 

4.0548 

1.2 

2.5482 

1.2296 

5.3579 

5.5099 

5.9728 

5.7169 

3.8902 

1.4 

1.6877 

.97768 

2.2061 

3.7225 

4.0550 

5.8902 

3.5813 

1.6 

1.1624 

.71458 

1.4461 

2.4855 

2.7265 

2.6295 

2.4319 

1.8 

.75837 

.46625 

.94741 

1.6407 

1.8157 

1.7669 

1.6481 

2.0 

.48740 

.30074 

.61629 

1.0813 

1.2125 

1.1948 

1.1269 

2.2 

.31426 

.19398 

.40500 

.72942 

.81971 

.81878 

. 78146 

2.4 

.19779 

.12233 

.26259 

.47637 

.55272 

.56177 

.54399 

2.6 

.12566 

.07802 

.17310 

.32108 

. 38025 

.39324 

.38603 

2.8 

.08288 

.05140 

.11710 

.22087 

.26599 

.27904 

.27716 

3.0 

.05376 

.03347 

.07919 

.15284 

.18783 

.20021 

.20133 

4.0 

.00467 

.00317 

.01135 

.02628 

.03541 

.04272 

.04585 

5.0 

.000923 

.000594 

.00258 

.00631 

.00921 

.01105 

.01207 

A * 

A 1.6,4* 

A * 

A 1.8,r* 

A » 

A 2.0,r* 

A * 

A 2.2,r* 

A 1 

A 2.4,r* 

A * 

A 2.6,r 

' A T 2 . 8 , r ’ 

A ♦ 

' 3.0,r 

1.1624 

.75852 

.48739 

.31425 

.19779 

.12561 

.08288 

.05376 

.71458 

.46625 

.30074 

.19398 

.13235 

.07802 

.05141 

.03547 

1.4461 

.94744 

.61628 

.40500 

.26259 

.17310 

.11710 

.079195 

2.4855 

1.6407 

1.0813 

.72042 

.47636 

.32109 

.22087 

.15284 

2.7265 

1.8156 

1.2125 

.81  971 

.55272 

. 58050 

.26601 

.18783 

2.6295 

1.7669 

1.1948 

.81878 

.5  6177 

.39325 

.27904 

.20021 

2.4319 

1.6481 

1.1269 

.78144 

.54399 

. 38608 

.27718 

.20133 

2.1759 

1.5196 

1.0485 

.75403 

.51634 

.37062 

.26848 

.19674 

1.5196 

1.5903 

.96677 

.68225 

.48479 

.35049 

.25572 

.18867 

1.0486 

.96679 

. S9099 

.63267 

.45270 

.32929 

.24190 

.17908 

.73403 

.68227 

.63266 

.58643 

.42132 

. 50816 

.22639 

.16382 

.51684 

.48479 

.45270 

.42182 

.39006 

.23592 

.21115 

.15751 

.37062 

.35050 

.32923 

.30816 

.28592 

. 26411 

.19612 

.14656 

.26847 

.25572 

.24153 

.22689 

.21115 

.19612 

.13072 

.13524 

.19674 

. 18866 

.17907 

.16382 

.15751 

.14656 

.15524 

.13427 

.04638 

.04580 

.04432 

.04240 

.03997 

.03628 

.03472 

.03202 

.01246 

.01245 

.01216 

.01169 

.01108 

.01041 

.00967 

.00893 
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Each  value  of  Afr,rf  listed  here  whould  be  divided  by  4tt 
since  the  functions  employed  had  been  normalized  to  1 instead 
of  1/4TT . To  obtain  the  values  of  Arr1  for  a particular 
energy,  the  multiplication  indicated  in  equation  (37)  was 
performed.  The  values  of  d)  (rT)  used  were  those  obtained 
as  mentioned  before  by  solving  the  equation  without  the 
exchange  integral.  Then  the  values  of  the  integral  at  the 
r chosen  is  obtained  by  graphing  the  values  of  the  corresponding 
Arr*  and  determining  the  area  mechanically.  For  example  the 
values  of  A^0,rT  for  = 2.1  yield  the  curve  drawn  in  Figure  3. 
It  should  be  noted  that  the  curve  was  not  carried  out  to 
rioc’but  rather  stopped  at  r-3.0  and  extrapolated  linearly  to 
r=6.0  whereas  the  actual  curve  would  extend  along  the  dotted 
lines.  This  definite  stop  and  linear  extrapolation  was  an 
attempt  to  correct  for  incorrect  distribution  of  the  Hartree 
4d  functions.  The  extrapolation  has  the  effect  of  pushing 
the  maximum  of  the  distribution  curve  for  the  4<lelectrons 
nearer  the  nucleus.  The  Hartree  functions  employed  here  did 
not  include  exchange,  i.e.  they  were  not  solutions  of  Fock 
equations.  Solutions  of  Fock  equations  for  Cu^ yielded  a 
shift  in  distribution  of  the  3d  electrons  in  the  same  direction 
as  the  linear  extrapolation  described  above.  A definite 
substantiation  of  the  writer* s method  of  integration.  Possibly 
this  accounts  for  the  fact  that  the  cohesive  energy  of  Ag  as 


30.  D.P.  Hartree,  cf.  ref.  21 
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computed  by  the  writer  compares  more  favorably  with  the 
observed  value  than  does  that  which  Fuchs  calculated  for  Cu. 

Now  to  obtain  the  cohesive  energy  one  must  assume  different 
values  of  the  energy,  solve  the  S - equation  with  exchange  in 
each  case,  apply  the  boundary  conditions,  and  determine  the 

corresponding  values  of  r0.  Back  in  Figure  3 the  curves  for 
<5=2.1  were  drawn.  The  dotted  lines  to  the  curve  with 
exchange  indicate  howr  the  boundary  conditions  are  applied 
graphically.  The  last  three  sentences  were  written  down 
easily  enough  but  actually  performing  what  they  required 
entails  solving  two  different  second  order  differential 
equations  by  numerical  integration,  two  hundred  and  twenty-five 
multiplications,  plotting  fifteen  curves,  evaluating  mechanically 
the  area  under  each  curve,  and  determining  the  values  of  r0 
by  plotting  the  curve  with  exchange  and  drawing  tangents  as 
described  above.  A.11  this  for  each  of  the  five  energies 
assumed.  It  is  true  in  this  discussion  that  it  has  been 
assumed  there  is  no  nuclear  motion  i.e.  absolute  zero.  However 
the  electrons  still  are  in  motion  and  the  energy  associated  with 
them  is  called  the  Fermi  Energy.  The  total  energy  of  the 
valence  is  the  sum  of  the  Fermi  energy  and  <S . The  Fermi 
energy  was  obtained  by  using  the  formula 
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Curve  I in  Figure  £ represents  the  energy  of  the  valence 
electron  without  Fermi  energy;  curve  II  the  energy  with  the 
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' 


29. 


Fermi  energy  added  and  the  dotted  line  the  ionization  energy. 
The  cohesive  energy  equals 


41. 


where  £ is  the  negative  of  the  ionization  energy  and  £ois  the 
energy  of  the  valence  electron  without  the  Fermi  energy.  This 
sum  can  be  obtained  directly  by  subtracting  the  energy 
corresponding  to  the  minimum  in  Curve  II  from  £ . The  lattice 
constant,  a^(edge  of  the  fundamental  cube)  is  related  to  rm,  the 
value  of  r corresponding  to  the  minimum  point  in  curve  Il^by  the 
equation 
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where  rm  is  in  atomic  units. 

The  values  of  the  cohesive  energy  and  lattice  constant  of 
Ag  as  determined  by  the  writer  are  given  in  Table  IX. 

Table  IX 

Calculated  Observed 

Cohesive  energy 

in  Kilo-cal/mol  95. Q 64.5 

Lattice  constant 

in  A°  2.72  4.07 

C . Conclusion: 

Since  the  calculated  results  for  the  cohesive  energy  and 
lattice  constant  for  Ag  are  distinctly  of  the  right  order  of 
magnitude,  the  purpose  of  the  dissertation  can  be  clearly  seen 
to  be  fulfilled.  While  the  lattice  constant  value  as  calculated 
by  the  writer  is  about  22%  worse  than  that  for  Cu  calculated 
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by  Fochs,  the  cohesive  energy  for  Ag  is  about  18$  better. 

In  a future  paper  the  writer  will  attack  the  problem  of 
further  refinement  of  this  method  in  an  endeavor  to  bring 
the  calculated  results  still  closer  to  the  observed  ones. 

The  reassurance  in  the  method  that  this  paper  has  given  will 
be  especially  gratifying. 
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Comprehensive  Digest 


The  credit  for  the  first  quantitative  attempt  to  account 
for  the  cohesive  energy  and  lattice  constant  of  solids  by  means 
of  wave  mechanics  belongs  to  Wigner  and  Seitz  whose  first  work 
appeared  in  1935.  They  worked  on  the  alkali  metal  Na  and  used 
the  semi-empirical  values  of  the  potential  of  Na  as  developed 
by  Prokofjew.  Their  results  agreed  well  with  the  experimental 
values  and  gave  definite  promise  for  continued  application  of 
the  Schroedinger  equation  to  solids. 

Accurate  potentials  of  the  type  employed  by  Wigner  and 
Seitz  for  elements  heavier  than  Na  are  not  available  and  recourse 
is  generally  made  to  the  results  of  Hartree  and  his  co-workers . 
The  use  of  Hartree  functions  necessitates  a form  of  the  radial 
wave  equation  for  the  valence  electron  which  is  much  more  com- 
plicated than  that  used  by  Wigner  and  Seitz.  The  reason  for  this 
is  that  the  exchange  interaction  between  the  electrons  is  absent 
from,  the  Hartree  functions  but  is  present  implicitly  in  the 
Prokofjew  potentials.  Pock  has  shown  that  this  exchange  inter- 
action between  the  valence  electron  and  the  core  electrons  can 
be  included  in  an  equation  of  the  form 


Fuchs  used  Hartre  functions  to  compute  the  cohesive  energy 
and  lattice  constant  of  Cu.  In  order  to  bring  his  results  into 
closer  agreement  with  experiment,  Fuchs  found  it  necessary  to 
appljy  several  correction  terms  to  his  original  determinations. 


1. 


(j(vr,v')  ^ znrir'e 


Because  of  the  addition  of  those  corrective  terms,  it  was  felt 
that  the  method  might  not  be  inherently  sound,  and  that  an  exten- 
sion to  heavier  metals  might  show  a complete  breakdown.  The  break- 
down would  manifest  itself  by  giving  completely  wrong  orders  of 
magnitude  for  the  cohesive  energy  and  lattice  constant.  Since  Ag 
is  a metal  closely  akin  to  Cu  in  structure  and  properties  but  has 
47  electrons  as  compared  with  the  29  of  Cu,  it  was  deemed  suit- 
able for  further  testing  the  validity  of  Fuchs T method.  Hartree 
functions  for  Ag-*-  have  been  calculated  by  Black  and  were  used  by 
the  writer  in  his  computations. 

Fuchs1  method  is  essentially  as  developed  in  the  succeeding 

paragraphs.  By  means  of  construe  ing  plane  at  the  midpoints  of 

the  lines  joining  an  atom  in  the  crystal  with  its  nearest  and 

next  nearest  neighbors,  a neutral  polyhedron  is  formed.  To  an 

approximation  that  has  be  n shown  to  be  well  within  the  tolerable 

limits  of  error,  the  polyhedron  may  be  replaced  by  a sphere  of 

equal  volume.  The  radius  of  this  snhere  is  called  ”r  ".  Then  since 

o 

each  sphere  is  nearly  neutral  one  can  concentrate  on  the  energy 
associated  with  one  and  neglect  the  rest. 

The  boundary  condition  for  the  wave  function  of  an  atom  in 
a crystal  has  been  shown  by  Bloch  to  be  that  the  wave  function 
must  be  smooth  as  it  passes  from  one  sphere  to  another.  This 

<Li 


-0  where ^ 


condition  id  expressed  mathematically  by  setting 

fy  ‘a, 

is  the  radial  wave  function  for  the  valence  electron.  If  one 
Ip  d.  (a>) 

sets  7 ' then  the  boundary  condition  is  transformed  to 

As 

This  condition  can  be  shown  to  be  satisfied  at  points  on  the  (fi 
vs.  r curve,  beyond  the  position  of  the  outer  shell  of  the  core 
electrons,  where  lines  drawn  from  the  origin  to  the  curve 


are  tangent.  Two  such  points  of  tangency  will  be  observed  at 
values  of  the  energy  numerically  less  than  that  of  the  stable 
energy  of  the  sphere.  At  this  latter  point  there  will  only  be 
one  tangent  and  at  points  corresponding  to  energies  greater  than 
this  stable  energy  no  tangent  can  be  drawn. 

Equation  (l)  is  the  Fock  radial  equation  for  the  valence 
electron  expressed  in  fundamental  units.  It  has  proved  more  con- 
venient in  this  particular  instance  to  employ  a different  set  of 
units  cabled  atomic  units.  If  these  units  are  used,  W set  , 

and  V set  -vc»the  equation  assumes  the  simpler  form 


equation  (2)  is  now  in  double  Rydberg  units.  This  equation  must 
be  solved  subject  to  the  boundary  condition  previously  discussed. 
The  standard  method  of  solution  in  such  cases  is  to  treat  the 
exchange  integral  as  a correction  term  and  neglect  it  in  the  first 
approximation.  Then  the  remainder  of  the  equation  is  solved 
by  numerical  integration.  The  value  of  (^thus  obtained  is  there- 
fore the  first  approximation  to  a solution.  In  the  alkali  metals 
the  value  thus  obtained  is  a fair  approximation  to  the  solution 
but  it  is  not  even  moderately  good  for  the  noble  metals.  The  rea- 
son is  that  the  interatomic  distance  is  much  larger  in  the  alkali 
metals  and  hence  the  exchange  effect  between  the  valence  and  core 
electrons  is  less.  Having  obtained  a value  for  (p  , we  may  now 
calculate  the  exchange  term.  As  it  stands  in  equations  (1)  and 


d r*  < ' -/.  ^ ' 

It  should  be  carefully  noted  that  the  energy  parameter, & , in 


it  does  not  lend  itself  particularly  well  to  numerical  calcu- 
lation. Fuchs  has  shwon,  however,  that  G(r,rT)  may  be  expressed 


in  another  form  as  in  equation  (5) 


5. 


where 


In  this  expression  the  functions  f(r)  and  f(rf)  are  the  Hartree 
radial  wave  functions  for  the  core  electrons  as  given  by  Black 
and  normalized  to  f^ir  . It  must  be  noted  that  the  summation 
extends  over  all  core  electrons  having  the  same  spin  as  the  valence 
electron,  that  is  over  one  half  the  number  of  core  electrons. 

Thus  the  exchange  integral  can  be  put  into  the  form 


energy  parameter.  A list  of  the  values  of  AT  . once  compiled 

P . P 
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For  clarity  in  discus  ion  it  has  proved  advantageous  to  set 


5. 


and 


6 . 


The  values  of  A1 ^ ^ f are  easily  seen  to  be  independent  of  the 


forms  a basic  table  from  which  the  values  of  A . , which  do 


T,r 


depend  upOh  the  energy,  may  be  computed  with  some  ease.  The 
exchange  integral  may  now  be  expressed  as 


7. 
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The  value  of  (7)  at  a -particular  r may  be  obtained  by  performing 
the  expansion  indicated  in  (5),  the  multiplication  indicated  in 
(6),  and  integrating  (7)  graphically.  Now 'having  calculated  the 
values  of  the  exchange  integral  at  several  values  of  r,  the  results 
can  be  inserted  in  (2)  and  the  now  inhomogeneous  equation  solved 
again  for  (j)  . This  0 is  then  plotted  against  r and  the  values  of 
rQ  corresponding  to  the  assumed  energy  plotted  on  an  <£,vs.  rQ 
graph.  Different  values  of  <£  are  chosen  and  the  process  repeated 
over  and  over  until  finally  an  £,is  assumed  at  which  there  is  only 
one  tangent.  This  value  of  C is  the  stable  energy  of  the  sphere. 

Thus  far  the  energy  of  motion  of  the  electrons,  the  Fermi 
energy,  has  been  neglected.  Since  the  purpose  of  this  paper  was 
to  test  the  validity  of  Fuchs*  method,  the  formula  that  Fuchs* 
used  for  the  Fermi  energy  will  be  employed  as  given  in  equation  (8) 


When  this  energy  is  added  to  the  6 vs.  r curve,  a new  minimum 

o 

is  obtained  at  rQ  - rm.  The  death  of  this  minimum  below  the  neg- 
ative of  the  ionization  energy  is  equal  to  the  cohesive  energy. 
The  value  of  r in  atomic  units  is  related  to  the.  a,  the  lattice 
constant  for  a face  centered  cubic  crystal,  by  the  equation 


The  cohesive  energy  and  lattice  constant  for  Ag  as  determined  by 
the  writer  are  given  in  Table  I. 


. 
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Table  I 

Calculated  Experimental 

Cohesive  energy- 

in  95.8  64.5 

kilo-cal/mol 

Lattice  constant 

in  2.72  4.07 

angstroms 

Since  the  results  are  distinctly  of  the  right  order  of 
magnitude  and  agree  with  the  experimental  values  to  the  expected 
extent,  Euchs*  method  can  safely  be  said  to  be  valid.  In  a future 
paper  the  writer  will  attack  the  problem  of  further  refinement 
of  this  method  in  an  endeavor  to  bring  the  calculated  results 
still  closer  to  the  observed  ones. 
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